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We discuss the relaxation of the scalar and pseudoscalar 
condensates after a rapid quench from an initial state with 
fluctuations. If we include not only the zero-mode but also 
higher modes of the condensates in the classical evolution, 
we observe parametric amplification of those "hard" modes. 
Thus, they couple nonlinearly to the "soft" modes. As a con- 
sequence, domains of coherent 7r-field emerge long after the 
initial spinodal decomposition. The momentum-space distri- 
bution of pions emerging from the decay of that momentum- 
dependent condensate is discussed. 



The exciting speculation that the dynamics of spon- 
taneous breaking of chiral symmetry could reorient the 
condensates in QCD as compared to the physical vac- 
uum has stimulated many studies, see e.g. and ||] 
for reviews. However, for the chiral symmetry breaking 
process in strong interaction the spinodal decomposition 
after a rapid "quench" does not yield a large correlation 
length II . Accordingly, a simulation Q of the classical 
equations of motion of the linear a-model on a 3d-lattice 
with spacing a = 1 fm confirmed that in the strong- 
coupling case no large domains form in which the 7r-field 
were essentially uniform, see also ||. 

We shall confirm that result but also argue that if the 
cutoff for the Fourier modes of the classical order parame- 
ter of chiral symmetry is pushed to larger values by some 
physical mechanism, e.g. parametric resonance [[To] [12j, 
domains much bigger than the Compton wavelength of 
the pion can form over time-scales 3> 1 fm. Thus, the 
"pseudoscalar laser" is a non-linear phenomenon devel- 
oping long after spinodal decomposition. 

Our rather simple idea is that parametric resonance 
leads to large occupation numbers of "hard" modes of the 
scalar and pseudoscalar condensates, which then couple 
non-linearly to the "soft" modes. Thus, the same mech- 
anism that is responsible for the explosive heating of the 
universe after inflation |l0| (see also were the inflaton 
decay was studied via a simulation on the lattice, similar 
to our studies here) leads to the formation of domains 
where the pseudoscalar condensate exhibits coherent os- 
cillations with big amplitude around its vacuum value. 
Those domains could be produced in the laboratory in 
collisions of protons |l|] or heavier nuclei Q at high ener- 
gies. 



As an effective theory of the chiral symmetry breaking 
dynamics J14| we apply the linear cr-modcl [fl5| 



+ ^(d^a + d^w) - U{a, t?) 



(0.1) 



The potential exhibiting the spontaneously broken sym- 
metry is 
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Here q is the constituent-quark field q — (u,d). The 
scalar field a and the pseudoscalar field ir = (^1,^2^3) 
together form a chiral field $ — (a, ir) . The parameters 
of the Lagrangian are usually chosen such that the chiral 
SUl(2) <8> SUr(2) symmetry is spontaneously broken in 
the vacuum and the expectation values of the conden- 
sates are (a) = f n and (7?) = 0, where f v = 93 MeV is 
the pion decay constant. The explicit symmetry break- 
ing term is determined by the PCAC relation which 
gives H = j^m^, where = 138 MeV is the pion 
mass. Then one finds v 2 = / 2 — to 2 /A 2 . The tr-mass, 
m 2 = 2A 2 / 2 + to 2 , which we set to 600 MeV, yields 
A 2 « 20. For the moment we focus exclusively on the 
dynamics of the chiral condensates, i.e. we set g = 0. We 
return to the case g 7^ at the end. 

We shal l di scuss the evolution of the system in the 
potential ( p.2| ), starting from a chirally symmetric ini- 
tial state, ($) = 0. Such an initial state could 
be reached thanks to the slow dynamics of the long- 
wavelength modes [^-^,^1; or, in the case g > 0, due 
to a possible first-order symmetry breaking phase tran- 
sition which "traps" the expectation value of the fields in 
the local minimum around ($) =0 i). We include fluc- 
tuations in the initial state by distributing the fields as 
Gaussian random variables with 



(0.3) 



($ 2 )=« 2 /16 , ($ 2 ) = v 2 /4 fm 2 



Since ($) = 0, the system is in an unstable condition 
for the potential at hand, eq. (0.2). The fields "roll 
down" towards the vacuum. In analogy to condensed 
matter physics this process is called spinodal decomposi- 
tion §,§. 
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FIG. 1. Time evolution of the 7r-field in the longitudinal 
space-direction, and isospin-0 direction, averaged over the 
transverse directions, time — 4 fm (practically independent 
of lattice spacing as long as 0.1 fm< a < 1 fm). 

We first show numerical solu tion s of the classical equa- 
tions of motion obtained from (0.1) on a lattice of physi- 
cal size 10 fmx 10 fm in "transverse" directions, and 40 fm 
in longitudinal direction. (Throughout the manuscript 
Minkowski- metric has been employed). We have inte- 
grated the equations using a second-order leap-frog al- 
gorithm, ensuring that energy and momentum were con- 
served to within 0.1%. We performed simulations with 
lattice spacing a = 1 fm as well as a — 0.2 fm. In the 
latter case, a coarse-graining procedure was applied to 
the initial configuration to obtain a correlation length of 
1 fm; this wa s done without artificially cooling the fluc- 
tuations (|5~5|). 

Fig. [l] depicts the 7r-field averaged over the transverse 
plane, showing the evolution within the first 4 fm. One 
observes that the field amplitude within each individ- 
ual small initial "domain" gets amplified, exactly as it 
should be ||. However, the wavelengths of the 7r-field 
oscillations more or less correspond to those of a free 
pion. Even for the rather small initial fluctuations ( |0.3| ) 
one can hardly observe big domains emerging during the 
initial "roll-down" HH, the less the larger the fluctua- 
tions at t = Q. Also, we have found that the amplitude 
of the 7r-field after "roll-down" are almost independent 
of a. This is no longer true for the long-time evolution 
to be discussed next. 

For large lattice spacing no big domains can be ob- 
served even at much later times, t = 80 — 100 fm, as 
shown in Fig. || The field is basically exhibiting only 
small and short-range oscillations. 



FIG. 2. As Fig. § but at time 80 - 100 fm (lattice spacing 
a — 1 fm). 

However, for a = 0.2 fm, one can observe long-range 
oscillations with a much bigger amplitude, clearly ex- 
hibiting the strong amplification of the long-wavelength 
modes and, correspondingly, big domains (Fig.||). More- 
over, this observation remains true even for larger initial 
fluctuations. 




FIG. 3. As Fig. but at time 80 - 100 fm (lattice spacing 
a = 0.2 fm). 

In Fig. |] we show the momentum-space distribution of 
7r-particles in the condensate (the distribution observed if 
the field decayed into particles at that moment) obtained 
from the spacial Fourier transform of ir(t,x). For sim- 
plicity, when computing dN/d 3 k we assumed the Hartree 
dispersion relation f2 2 = k 2 + A 2 ((er 2 ) + 5(7r 2 )/3 — v 2 ). 
(The even more simple-minded application of the free- 
pion dispersion relation gave practically the same result.) 
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FIG. 4. Distribution of pions of given isospin in mo- 
mentum space, dN/d 3 k (lattice spacing a — 0.2 fm, time 
t = 100 fm). The broken line depicts a statistical (Bose) 
distribution. 

The result has been averaged over several initial field 
configurations, chosen randomly according to (3.3). One 
observes a very sharp rise below k ~ 70 MeV, proving 
the spacial homogeneity of the pseudoscalar condensate 
at t = 100 fm. Note that a statistical (Bose-Einstein) 
distribution does not show such a "rise" , i.e. local max- 
imum of the second derivative, at non-vanishing k. (In 
Fig. |I| the statistical distribution is for a temperature 
T = 160 MeV and a chemical potential fi n = 100 MeV.) 

To understand the results of the simulation we analyze 
the effect of the cr-field oscillations on the Fourier modes 
of the 7r-field in a more schematic model. We employ the 
linearized classical equation of motion for the pion field, 
neglecting also explicit symmetry breaking (m, = 0). 
Furthermore, we assume that the sigma-field makes har- 
monic oscillations with amplitude a$ in the minimum of 
the potential, a = aosm(m a t). After Fourier transfor- 
mation, the evolution of the 7r-field modes becomes (see 
also 0) 



7Tfc + {k 2 + A 2 (T (cr + 2v)} TT k = 



(0.4) 



The expression in curly brackets gives the frequency f2| of 



the 7r-field oscillations, 
we assume TTk(t = 0) 



As initial condition for the 7r-ficld 
= 1/VXh, = 0) = yfthfi. 



The solution of this equation can be found numerically 
without great difficulties (we used the routine "odeint" 
from modified to work with double precision, and 
required accuracy 10 -9 ). 

Fig. shows the occupation number 
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of the 7r-field modes for three different cr-field amplitudes. 
(For (To « f w found much smaller occupation num- 
bers and "squeezed" resonance bands.) Oscillations with 
amplitude <Jq > v will certainly not be harmonic, see 



e.g. p^ ]. However, this is only a schematic and quali- 
tative example anyway, since it completely neglects the 
nonlinearities in the 7r-field and the back-reaction on the 
cr, which are most essential in the lattice simulation [jl3| . 
We thus do not attempt to be more "realistic" in this 
respect. 

The strong amplification of the soft modes of the ir- 
field, k ~ has been interpreted as being due to an 
imaginary effective mass of the pion . However, there 
are additional resonance bands at higher k, being partic- 
ularly pronounced if the amplitude of the cr-field is big. 
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FIG. 5. Occupation numbers of various Fourier modes of 
the 7r-field after 10 oscillations of the cr-field with amplitude 
Co = v (thick line with very strong soft modes), <ro = 2u 
(thinner line), and ao ~ rn a (thinnest line), respectively. 

In Fig. U we show the energy density of the 7r-field, 
d 3 k 



e w = 3 



(2tt) 3 
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The factor 3 in front of the integral comes from the three- 
fold isospin degeneracy of the 7r-field. This is the energy 
density that would be measured if the classical field de- 
cayed into particles with energy Qk at the corresponding 
time. One observes that the energy density of the 7r-field 
increases exponentially in time. In a more realistic cal- 
culation, as e.g. the lattice simulation described above, 
that exponential growth will shut off eventually due to 
energy conservation and the back-reaction. 

In a case like this we do not integrate out the k > 
modes, or the modes above some cutoff k c <C T, and 
treat only the zero-mode classically |Q (evolving in the 
effective potential generated by the "hard" modes above 
k c ). In particular, we can not compute the dissipative 
corrections to the classical dynamics of the soft modes, 
which is due to scattering with the "hard" modes, per- 
turbartively (for an alternative approach see H|). We 
therefore decided to describe all modes within the first 
few parametric resonance bands classically and to solve 
the full nonlinear field equations (on the lattice). 



Of course, we can not include arbitrarily hard modes 
in our classical description as we would run into the well- 
known Rayleigh- Jeans type of ultraviolet divergences, cf. 
e.g. the discussion in Fortunately, our lattice simu- 
lations were practically independent of the lattice spac- 
ing (domains always formed) as long as a ~ 0.1 — 0.3 fm, 
showing that higher resonance bands do not contribute to 
the classical evolution. One may have to employ smaller 
a (i.e. larger k c ) i f the initial fluctuations were substan- 
tially larger than ( |0 . 3| ) - 
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Furthermore, we found that the cr-field relaxes to a 
value more closely to /„•, the correct vacuum value, if a 
is small. We obtained (a) = 0.80/^ for a = 0.2 fm while 
(<t) = 0.56/tt for a = 1 fm; both values determined at 
t = 120 fm; see Fig. ||. This observation is opposite to the 
expectation in the static, thermal equilibrium case |L7| . 
In the dynamical case the high-frequency classical fluctu- 
ations, which would otherwise mimic finite-temperature 
contributions to the potential, cool more efficiently by 
transfering more of their energy to the soft modes if a 
is small (as is obvious from Fig. || and ||). This is not 
in contradiction to the fact that for t — > oo one should 
have equipartitioned energy (classically) as the reverse 
process, energy transfer from soft to hard modes, occurs 
on much longer time-scales. Thus, we employ k c ~ 1.5m CT 
in the simulation of the model on the lattice (described 
above). In this way, we account for the full non- linear 
evolution of those parametrically amplified modes. 

The modes above k c could be accounted for by inte- 
grating them out of the partition function, assuming e.g. 
a thermal distribution |^,[I7|. The scattering of modes 
with k < k c off modes with k > k c will then introduce 
fluctuation and dissipation terms into the classical equa- 
tions of motion. Alternatively, the interactions between 
classical modes and quantum fluctuations could even be 
treated explicitly, such that the latter effectively act as a 



heat-bath for the former [ pj| |. That is, however, beyond 
the scope of the present manuscript. 

We now reintroduce the coupling of the cr-fiel d to the 
quarks. For the constituent quark mass m q = gV& 2 + tt 2 
in the vacuum to be one third of the proton mass, we 
employ g = 3.3. Furthermore, dropping the 7r-field and 
again assuming a homogeneous cr-field, one arrives at the 
following equation of motion for modes of the auxilliary 
field V 



i>k + {k 2 + g 2 [o + v f - ig&} ip k = 



(0.7) 



where a = &o sinm CT i, as above. The quark wavefunction 
can be obtained from tpk as 



1% t) = [*7 ' d + 9°\ J 
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R r denote the two eigenvectors of 7 with eigenvalue +1. 
Aga in, t he real part of the expression in curly brackets in 



eq. (0.7) denotes the oscillator frequency £l 2 . As initial 



condition we employ 
^ k (t = Q)=Af k 



Mt = 0) = -mkMk 



(0.9) 



with 1/A4 = y2f2fe (flk + g& + gv). The occupation 
number of each mode is 
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2n k 
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and the energy density is given by expression (0.6), with 
a factor 24 (instead of 3) in front of the integral. 
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FIG. 7. Occupation numbers of various Fourier modes of 
the i/j-field after 10 oscillations of the cr-field with amplitude 
— v (bottom curve) and ao = m a (top curve; 1 + rJ k is 
shown in that case). 
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Fig. |7j shows the occupation number of the various 
fcrmion modes. Due to the Pauli principle exponential 
growth of the occupation numbers is not possible. How- 
ever, while the incoherent process a — > qq with "on-shell" 
kinematics is forbidden (the quark mass, averaged over 
one oscillation of the er-field, is m q = g(f n + 2ao/ir) > 
m a /2), the cr-modes can act as coherent field to produce 
fcrmions. This is because the quark-mass becomes small 
during some time-interval of the coherent cr-oscillations 
and the frequency flk of the quark wavefunction can drop 
below m<T (one can show analytically that no quarks are 
being produced if the tr-field is static relative to Vfe)- 
As opposed to the bosonic occupation numbers, small-fc 
quarks are produced if the amplitude of the driving a- 
field is large. However, even in this very simple analysis 
disrespecting total energy conservation the long-time en- 
ergy density of the quarks is well below that of the pions. 
Only if the cr-oscillations were strongly damped, e.g. by 
a rapidly expanding background, quark-pair production 
may be as important as pion production. The discussion 
of the condensate relaxation in an expanding metric is 
out of the scope of the present manuscript, see e.g. HQ. 

In summary, we have performed numerical simulations 
of the classical dynamical evolution of the chiral conden- 
sates in strong interaction on a lattice. In particular, we 
have focused on the situation where the chiral fields are 
rapidly quenched from an initial state with fluctuations. 

We confirmed previous results in that the spinodal de- 
composition process does not lead to formation of big 
domains of coherent 7r-ficld. This remains true also at 
later times if indeed only the zero-mode can be described 
classically, i.e. for rather large lattice spacing a > 1 fm. 
However, if higher modes are also treated classically by 
decreasing the lattice spacing to ~ 0.1 — 0.2 fm, long- 
range correlations (domains) do form as a consequence 
of parametric resonance and the nonlinear dynamics of 
those exponentially amplified modes. Despite evident 
limitations of the effective model that we employed, we 
think that the strong amplification of "hard" modes of 
the condensates, which as a consequence evolve nonlin- 
early and eventually drive the amplification of the "soft" 
modes, may occur also in more realistic models. 

Finally, we have studied very briefly the effect of cou- 
pling the chiral fields to constituent quarks. Despite the 
fact that the Dirac vacuum is not stable in the presence of 
coherent classical a-field oscillations with big amplitude 
or high frequency, the energy density of the quarks re- 
mains well below that of the 7r-field modes at late time. 
Thus, the coupling of the condensates to quarks is not 
likely to destroy the "pseudoscalar laser" emerging if an 
initial state with fluctuations is rapidly "quenched" to 
the spontaneously broken chiral symmetry. 
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FIG. 8. Additional Figure: cr-field (averaged over the entire 
lattice) as a function of time for a = 0.2 fm (full line) and 
a — 1 fm (dotted line). 



